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ABSTRACT 


The  ability  to  predict  surface  currents  can  have  a  beneficial  impact  in  several 
activities,  such  as  Search  and  Rescue  and  Oil  Spill  Response,  as  well  as  others  more 
purely  scientific,  operational  or  economic.  The  Naval  Postgraduate  School,  in 
conjunction  with  the  Romberg  Tiburon  Center  and  the  University  of  California  Santa 
Barbara,  has  been  studying  this  purpose  for  the  Coastal  Response  Research  Center, 
University  of  New  Hampshire  and  NOAA.  So  far  the  prediction  was  based  on  tide  and 
persistence  of  the  reminiscent  current.  Faced  with  increasing  error  under  changing 
environmental  conditions,  further  study  of  other  influences  became  fundamental,  in  order 
to  increase  reliability.  This  study  is  a  part  of  that  effort  by  studying  the  impact  of  the 
wind-induced  currents  on  forecasting. 

Based  on  a  year  and  a  half  of  wind  and  HF  surface  current  readings,  the  wind 
surface  current  interaction  is  analyzed  and  quantified.  Then  that  influence  is  plugged  into 
the  forecast  algorithm.  The  final  results  show  that  the  wind  driven  surface  current  is 
about  2%  of  the  wind  magnitude  rotating  around  50°  clockwise,  with  coherence  after  17h. 
The  wind  introduction  into  the  forecast  improved  the  accuracy,  but  only  by  an  average  of 
10%.  The  error  still  climbs  with  the  variability  of  the  environment,  but  knowing  the  wind 
influence  allows  other  factors’  influences  to  be  observed  more  accurately,  such  as  the 
magnitude  of  the  current  itself.  Forecasting  is  now  done  with  0.15  m/s  plus  or  minus  0.1 
m/s  at  95%  confidence. 
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I.  THE  SURFACE  CURRENTS  FORECAST 


A.  APPLICATIONS 

The  use  of  the  sea  is  not  natural  to  man.  Instead  of  moving  with  a  solid  surface 
underneath  his  feet  and  an  almost  intangible  fluid  surrounding  him,  at  sea  he  has  to  deal 
with  a  much  more  viscous  fluid.  One  of  the  major  differences  that  this  produces  is  that  on 
earth,  if  he  stops  moving  his  body,  or  his  car,  he  stops.  At  sea,  he  doe  not.  In  order  to  take 
the  best  advantage  of  this  distinction,  understanding  how  this  movement  works  is 
necessary.  In  the  fifteenth  and  sixteenth  centuries,  the  Portuguese  voyages  of  discovery 
that  navigated  the  entire  globe  in  simple  sailing  ships  relayed  much  of  this  knowledge. 
Charts  of  currents  were  drawn  and  pilots  were  written  in  each  trip. 

With  the  passage  of  four  centuries,  several  approaches  have  evolved  to  face  the 
challenges  of  today.  One  of  them  is  the  remote  evaluation  of  the  surface  current  by  HF 
RADAR  (explained  more  in  detail  on  II.B.2.a.).  This  method  is  excellent  in  filling  the 
gap  between  satellite  remote  sensing  measurements  and  the  shore;  however,  it  is  not  good 
for  the  Surf  zone,  for  which  other  methods  can  be  applied,  like  the  one  described  in  Puleo 
et  ah,  2002.  Depending  on  the  frequency,  HF  radar  can  be  used  to  read  from  around  3  to 
200  km  off  the  coast.  Surface  currents  are  obtained  with  regular  readings  through  this 
means,  allowing  the  creation  of  a  real  time  surface  picture.  Several  works  are  still 
debating  the  accuracy  of  this  system  and  what  influences  it.  This  study,  focused  on  the 
forecast  of  the  next  24  hours  (24h)  surface  currents,  takes  the  next  logical  step.  Extracting 
currents  using  electromagnetic  radiation  (EMR)  on  the  HF  band  implies  that  all  the 
readings  made  are  of  the  first  two  meters  (2m)  of  depth,  because  EMR  does  not  reach  any 
further  down  and  in  this  shallow  layer,  wind  is  considered  to  be  a  major  driver  of  the 
currents. 

There  are  not  many  studies  considering  the  interaction  between  the  previously 
mentioned  surface  layer  currents  and  wind,  being  almost  all  based  on  short  time  series 
and  with  limited  resolution,  most  under  the  inertial  period  (24h).  The  existence  of  a 
database  of  two  time  series  of  both  wind  and  current,  with  lengths  over  one  year  and  with 
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an  hourly  resolution,  presents  an  opportunity  to  unveil  some  properties  of  that  interaction 
that  are  currently  poorly  known.  Based  on  these  time  series  data  (described  in  II. B.),  this 
thesis’s  aim  is  to  first  find  the  surface  current  component  that  is  driven  by  the  wind,  by 
constructing  a  transfer  function  between  them  (described  in  Section  II).  This  transfer 
function  will  then  be  used  as  an  extra  input  for  an  already  existing  forecast  model  that 
relies  solely  on  persistence  and  tides  (described  in  I.B.),  expecting  to  bring  better 
accuracy  and  reliability  to  it.  The  results  obtained  for  the  wind/surface  current  interaction 
in  this  long  time  series  may  also  be  applied  to  many  other  future  studies. 

Due  to  its  practical  implications,  several  institutions  are  interested  in  systems  that 
allow  this  knowledge  of  the  surface  currents  and  in  its  forecast,  because  it  has  direct 
impact  to  every  sea  user  that  needs  to  plan  any  type  of  short  term  deployment.  That  is  the 
case  of  the  National  Oceanic  and  Atmospheric  Administration  (NOAA)  Coastal 
Response  Research  Center,  who  funded  a  study  on  “Delivery  and  Quality  Assurance  of 
Short-Term  Trajectory  Forecasts  from  HF  Radar  Observations,”  of  which  this  thesis  is 
also  a  part.  Better  knowledge  of  the  next  24h  surface  currents  may  be  used  to  get  better 
results  in: 

1,  Navy  Type  Operations 

Surface  current  forecasting  has  a  significant  impact  on  non-motorized  activities, 
because  propulsion  reduces  the  dependency  on  the  flow  to  which  the  user  is  exposed. 

a.  Special  Operations 

Deployment  of  divers  or  any  type  of  raft  maneuvering  at  the  coast  may 
benefit  from  a  shortened  mission  time  and  energy,  reducing  risks  and  costs. 

b.  Mine  Warfare 

Locating  the  positions  of  drifting  mines  remains  a  modern  reality  because 
mines  are  still  the  cheapest  and  most  effective  way  of  making  constraints  to  the  use  of  a 
certain  sea  area  or  a  harbor.  To  avoid  or  to  sweep  drifting  mines,  knowing  where  a  field 
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is  today  will  not  help  much  tomorrow,  because  the  mines  drift  with  the  current  at  their 
depth.  If  you  are  deploying,  you  need  to  know  if  the  mine  behavior  will  be  what  you 
intend. 


c.  Flag  Presence 

Most  flag  presence  missions  consist  of  long  periods  where  a  ship  is 
constrained  to  a  certain  position  or  box.  Taking  advantage  of  the  currents  will  result  in 
better  efficiency  in  the  use  of  time  and  fuel. 

d.  Maneuvers  and  Evolutions 

To  maneuver  a  ship  a  long  distance,  the  vector  of  her  self  propelled 
velocity  must  be  added  to  the  surface  current  velocity,  also  resulting  in  better  efficiency 
in  the  use  of  time  and  fuel. 

e.  Weapons  Practice 

Target  drift  in  time. 

In  live  surface  fire  practice,  usually  a  target  is  deployed  or  towed  at  a  very 
slow  speed.  More  accurate  knowledge  of  the  current  will  predict  target  drift  over  time, 
allowing  a  better  residence  time  in  the  exercise  area  with  fewer  maneuvers,  making  it 
safer. 


2.  Coast  Guard  Type  Operations 
a.  Search  and  Rescue 

Most  distressed  platforms  do  not  have,  or  only  have  partial,  locomotion 
capability.  Examples  are  man  overboard,  life  saving  rafts,  debris  or  a  damaged  ship. 
Having  the  correct  surface  current  scenario  will  allow  two  major  advantages: 

(1)  Reduces  Position  Uncertainty  to  the  Distressed  Position. 
The  future  target  position  is  a  function  of  time  multiplied  by  the  elements’  influence. 
When  the  rescue  craft  reaches  the  search  zone,  it  usually  has  to  face  an  area  which  size 
depends  on  the  uncertainty  of  the  surface  current  and  the  wind  influence  on  the  target 
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(and  the  initial  position  error).  Reducing  surface  current  uncertainty  will  directly  reduce 
the  search  area,  reducing  the  time  needed  to  find  the  target,  allowing  better  survivability 
odds. 

(2)  Rescue  Missions  Savings  in  Hours  and  Fuel.  Added  to  the 
advantage  in  (1)  of  facing  a  shorter  area,  rescuers  are  able  to  plan  a  more  accurate 
intercept  course,  arriving  earlier  than  possible  previously  which  results  in  an  even  tighter 
search  area.  On  the  scene,  better  sweeping  plans  can  be  made. 

b.  Oil  Spill  Response 

Knowledge  of  the  impact  area,  where  a  spill  occurred  and  where  it  is 
going,  allows  the  efficient  deployment  of  the  available  means,  reducing  damage.  For  a 
good  result,  it  must  be  coupled  with  the  pollutant’s  dispersive  characteristics.  Also,  it 
helps  track  pollutants  that  may  have  come  from  tank  washings. 

c.  Illegal  Traffic  and  Immigration 

Much  of  this  traffic  is  done  by  small  slow  boats,  for  which  the  surface 
current  is  a  determinant  of  the  true  course.  Also,  sometimes  packages  are  sent  adrift  on 
purpose  or  to  avoid  capture.  Knowing  the  currents  will  facilitate  the  collection  of  them. 

d.  Patrol 

Better  efficiency  in  the  use  of  time,  fuel  and  area  covered. 

3.  Fisheries 

a.  For  Those  Who  Monitor 

One  typical  periodic  problem  is  the  understanding  of  how  a  bloom  will 
spread  out.  Knowing  more  about  buoyant  behaviors  helps  avoid  the  capture  and 
consumption  of  organisms  that  may  have  toxic  elements.  Also,  surface  current 
knowledge  gives  a  better  understanding  of  the  probable  paths  of  buoyant  (or  slow 
velocity)  organisms,  such  as  eggs,  larvae  or  youngsters. 
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b. 


For  the  Fisherman 


Coastal  fishermen  depend  a  lot  on  timing  and  their  profits  depend  a  lot  on 
fuel  consumption.  The  knowledge  of  how  the  currents  will  behave  provides  a  tool  to 
make  transits  and  staying  with  fisheries’  spots  more  cost  efficient.  They  also  benefit  from 
more  accurate  knowledge  of  the  species  behavior  described  in  the  previous  paragraph. 

4.  Meteorological/Oceanographic 

Other  than  the  direct  relevance,  it  is  an  important  tool  to  models,  as  they  need 
periodic  assimilation  of  reality  data  to  be  accurate.  This  means  can  also  control  model 
forecasts’  quality  by  direct  comparison  with  its  own. 

5.  Commercial  and  Recreational  Navy 

Coastal  transits  can  be  planned  to  be  more  efficient  in  terms  of  fuel,  time  and 

safety. 


6.  Habitat  Impact 

Knowledge  of  buoyant  behaviors  makes  it  possible  to  understand  the  behavior  of 
an  atypical  discharge  from  a  river  or  sewer  system,  providing  a  tool  to  better  manage  the 
situation. 

B.  ACTUAL  MODELS 

Surface  currents  are  usually  collected  by  direct  measurements  using  current 
meters  or  current  profilers,  or  by  remote  sensing  systems,  such  as  coastal  radars  or 
satellites.  In  order  to  make  forecasts,  there  is  need  of  actual  observations  plus  a  model  for 
future  behavior,  and  forecasting  ability  will  be  in  direct  relation  to  the  type  of 
observations  available. 

Retrieved  data  from  satellite  remote  sensing  is  usually  based  on  readings  of  sea 
surface  temperature  (SST),  or  other  pixel  differentiation  infonnation,  to  allow  tracking  of 
characteristics  by  maximum  cross  correlation.  This  gives  resolutions  good  enough  for 
coastal  applications  that  range  from  1  km  to  25  km  according  to  the  frequency  used,  the 

latter  for  microwave,  which  is  the  only  method  able  to  penetrate  clouds  (Emery  et  ah, 
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2005).  They  are  then  coupled  with  satellite  altimetry  infonnation  to  create  the 
geostrophic  flow  and  make  corrections,  if  a  satellite  is  already  on  task.  However,  there 
are  two  primary  disadvantages  to  using  satellites:  they  are  more  expensive  to  operate  and 
have  less  time  resolution,  despite  the  possibility  of  coupling  with  models  in  between. 
Since  the  time  resolution  is  twice  a  day,  a  function  of  the  orbit,  it  does  not  allow 
handling  of  smaller  time  scale  phenomena.  Those  phenomena  may  not  be  very  important 
in  the  open  sea,  but  they  gain  substantial,  if  not  major,  influence  near  shore.  One  of  such 
cases  is  the  tides.  Also,  tides  cannot  be  forecasted  with  poorer  resolutions.  No 
quantitative  evaluation  of  the  accuracy  of  these  models  was  found  in  the  references. 

The  use  of  radar  systems  for  this  purpose  is  becoming  widely  used  throughout  the 
world  because  of  the  many  unique  applications  that  exist  for  real-time  surface  current 
mapping  data.  Also,  it  is  a  cheaper  and  less  dependent  system  compared  with  the  use  of 
satellites,  particularly  if  the  user  is  interested  primarily  in  coastal  currents.  Since  the  real 
time  data  is  already  collected,  several  academic,  operational  and  commercial  institutions 
have  been  trying  to  go  further  and  forecast  the  surface  currents.  Following  are 
descriptions  of  two  models  referencing  this  work. 

1.  The  Naval  Postgraduate  School  Model 

This  approach  relies  on  the  Coastal  Ocean  Dynamics  Applications  Radar 
(CODAR  -  a  HF  RADAR  system)  readings  of  surface  currents,  as  described  in  II.  B.  2. 
a..  These  readings  are  then  treated  for  Geometric  Dilution  of  Precision  (GDOP),  which 
takes  into  account  the  uncertainty  provoked  by  the  angle  of  the  azimuths  between  stations 
and  reading  position  (the  closer  to  90  degrees,  the  better).  The  currents  associated  with 
tides  are  then  removed  by  software  from  MathWorks  called  “t_tide.”  More  information 
about  “t  tide”  can  be  found  in  Pawlowicz  et  ah,  2002,  which  performs  a  best  fitting 
frequency  relation  between  the  time  series  and  the  tidal  harmonic  constituents,  taking  the 
period  of  interest  into  account.  This  model  uses  a  Cook,  2008,  alteration  to  the  software. 
After  removing  the  tidal  effects  from  the  current,  the  remaining  current  is  then  averaged 
to  the  mean  current  for  that  day.  Based  on  persistency,  this  current  is  then  projected  to  the 
next  24  hours  and  is  added  to  the  tidal  current  that  is  again  forecasted  by  the  “t  tide” 
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prediction  software.  The  accuracy  of  this  model  is  typically  7km  distance  after  24h,  but 
frequently  exceeds  20km.(  Garfield  et  ah,  2007).  This  thesis  builds  on  this  model 
improvement  by  introduction  of  the  wind  influence. 

2.  The  Coast  Guard  /  University  of  Connecticut  Model 

Working  with  the  U.S.  Coast  Guard,  the  University  of  Connecticut  developed  the 
Short  Term  Predictive  System  (STPS,  Ullman  et  ah,  2003;  O’Donnell  et  ah,  2005)  for  the 
U.S.  Coast  Guard  Research  &  Development  Center.  This  model  also  relies  on  CODAR 
readings  treated  for  GDOP.  In  their  work  several  essays  were  attempted.  In  the  last,  the 
surface  current  was  divided  into  three  parts:  tidal  flow,  wind  driven  and  mean.  The  tides 
and  the  wind  influence  are  estimated  by  a  different  algorithm,  the  Gauss-Markov 
estimator,  which  divides  what  part  of  a  signal  is  tidal  or  wind  motivated  and  what  is  not, 
by  comparison  of  time  series,  creating  weights  that  are  calculated  from  the 
autocovariance  functions.  This  signal  decomposition  then  allows  prediction  of  the  next 
25h  by  persistence,  tide  and  predicted  winds  influence.  One  NDBC  buoy  is  used  for  the 
wind  data,  with  results  more  accurate  than  those  from  the  previous  Coast  Guard  results 
using  only  NOAA’s  tidal  data.  Root  Mean  Square  (RMS)  errors  were  in  the  order  of  10 
cm/s  (ten  centimeters  per  second)  for  each  component  of  the  velocity,  u  and  v 
(corresponding  to  14.14  cm/s  for  the  vector).  Forecasts  for  search  boxes  were  then  made 
with  20  cm/s  RMS  in  order  to  cover  uncertainty  in  all  measurements  and  then  compared 
with  real  data  from  drifters  and  statistical  approaches  by  a  Monte  Carlo  simulation. 

The  main  differences  between  the  Coast  Guard  model  and  the  one  used  in  this 
thesis  are:  the  length  of  the  evaluation  data  set  -  seven  days  versus  1,5  years,  the  way  the 
estimation  for  the  tides  and  winds  influences  are  calculated  -  Gauss  Markov  during  one 
Month  for  the  tides  and  during  1  day  for  the  wind  versus  t  tide  for  one  Month  and  time 
series  analysis  for  1 ,5  years  respectively,  the  wind  data  position  (number  of  buoys)  -  one 
versus  five  different  positions,  and  the  different  number  of  CODAR  data  stations  -  two 
versus  10  different  stations.  The  STPS  concludes  that  “...another  analysis  approach  to 
quantifying  the  wind-driven  circulation  is  called  for.” 
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C.  FORECAST  AND  THE  INTRODUCTION  OF  WIND  INFLUENCE 

All  the  works  concerning  forecast  of  surface  currents  point  to  the  necessity  of 
accurately  knowing  all  its  dependencies  and  their  magnitudes.  We  know  that  the  end 
result  will  be  in  direct  proportion  to  the  accuracy  of  the  data  and  of  the  relations  we  plug 
in.  One  of  those  dependencies  is  the  wind,  as  Ekman’s  work  proved  by  taking  one  step 
further  our  intuitive  experience  that  the  wind  influences  the  movement  of  the  water 
underneath  it  (Ekman,  1905). 

The  work  reported  here  was  done  after  the  initial  results  of  the  use  of  the  Naval 
Postgraduate  School  model  for  the  NOAA’s  Coastal  Response  Research  Center  (Garfield 
et  al,  2007),  where  the  need  to  understand  the  wind’s  influence  on  the  surface  current  was 
identified  as  possibly  relevant,  among  other  factors.  Despite  this  thesis’s  focus  on  the 
wind,  other  dependencies  will  also  be  analyzed  when  doing  error  assessment.  In  order  to 
reach  a  forecast  model  including  the  wind,  a  full  time-series  analysis  study  of  wind  / 
surface  current  interaction  is  made  and  shown  in  section  II.  After  detailing  this 
interaction,  and  having  a  transfer  function  created,  a  forecast  model  is  put  to  work 
producing  results,  which  will  be  described  in  Section  III.  Conclusions  from  these  results 
are  then  extrapolated  and  discussed  in  Section  IV. 
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II.  THE  WIND  /  SURFACE  CURRENT  INTERACTION 


A.  PREVIOUS  WORKS  USED  IN  THIS  STUDY 

Despite  the  long  list  of  studies  on  air  /  sea  interaction,  few  studies  were  found 
concerning  strictly  the  interaction  between  the  wind  and  the  surface  current,  with 
quantified  results  for  periods  less  than  24  hours.  Nevertheless,  portions  of  the  following 
works  contribute  to  the  understanding  of  the  wind  surface  current  relation  process: 

1.  Ekman,  1905 

Although  this  thesis  is  based  on  mainly  empirical  results,  Ekman  provided  the 
basic  theoretical  support  for  it  (Ekman,  1905).  Starting  with  the  geostrophic  equations 
with  vertical  friction  and  dividing  the  flow  into  geostrophic  and  frictional  components,  he 
related  the  wind  stress  and  its  influence  on  the  oceanic  flows  as  a  function  of  depth, 
creating,  among  other  results,  what  we  now  know  as  “Ekman  Spiral”  (Figure  1).  The 
Ekman  theory  assumes  infinite  ocean  and  constant  eddy  viscosity  (Av),  which  is  an 
approximation  of  an  already  poorly  known  value  (Radko,  NPS  class  notes,  2007). 


Figure  1.  Ekman  Spiral.  -  The  angle  between  the  wind  and  the  surface  current  is  45° 
clockwise  in  the  Northern  Hemisphere  (From:  http://oceanservice.noaa.gov, 
September  2008). 
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According  to  his  work,  the  wind  stress  is  given  by: 


.  dv 

r  =  p0A  — 

dz 

where  r  is  the  wind  stress  and  p0  represents  the  water  density  as  a  constant. 


The  amplitude  in  depth  is 


\Ve  1=1  vo\e 


f  being  the  Coriolis  term  and  z  the  depth. 
The  direction  is  shown  by: 
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cp{)  being  the  angle  at  the  surface. 


(1) 


(2) 


(3) 


For  this  study,  the  two  main  results  are: 


•  The  angle  of  the  current  at  the  surface,  which  is  45°  to  the  right  of  the 
wind  stress. 

•  The  amplitude  of  the  velocity  shown  by 
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2.  McNally,  Luther  and  White,  1988 

This  is  a  study  on  subinertial  frequency  response  of  wind  driven  currents  in  the 
mixed  layer.  The  data  set  consisted  of  63  tracks  of  drifters,  between  1976  and  1981,  with 
results  being  binned  at  5  days  periods.  The  drifters  operated  at  a  30m  depth.  The  study 
described  a  process  of  the  interaction  between  wind  and  currents,  despite  the  different 
objective,  analyzed  depth  and  periods.  Among  the  results  are  descriptions  of  the 
following  effects  indirectly  related  to  this  work. 
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a.  The  angle 

The  angle  was  found  to  be  around  25°  to  the  right  of  the  wind  stress  (at 
30m),  which  means  the  existence  of  the  clockwise  rotation  in  depth,  despite  the  angle, 
was  smaller  than  the  expected  by  Ekman  theory. 

b.  Correlation 

The  correlation  was  found  to  be  significant  between  near  inertial  to  16 
days.  No  relation  was  found  under  the  inertial  period. 

c.  Major  Axis 

This  was  found  to  be  from  75°  to  the  right  of  the  wind  stress  at  near 
inertial  to  15  0  at  the  lower  frequencies. 

d.  Lags 

Lags  were  observed  from  30  0  at  near  inertial  to  zero,  corresponding  to 
under  four  hours  in  time. 

3.  Foster,  1993 

This  thesis  studied  the  relation  between  diurnal  surface  winds  and  currents  in 
Monterey  Bay.  The  author  analyzed  the  daily  wind  surface  relation  for  September  1992 
in  Monterey  Bay,  using  two  CODAR  sites  for  the  currents,  two  moorings  and  three 
coastal  stations  for  the  winds.  The  result  was  data  on  the  behavior  of  the  wind  and  current 
circulation  and  interaction  in  this  zone  of  the  Monterey  Bay.  The  major  results  that 
impact  this  study  are  the  current  response  to  the  breezes  be  under  two  hours  (the  time 
resolution  of  the  data)  and  the  fact  that  the  more  offshore  currents  rotate  clockwise  in 
time,  to  the  right  of  the  wind. 
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4.  Xu  and  Bowen,  1993;  Monismith  et  al.,  2006 

This  is  a  work  on  wave  and  wind  driven  flow  interaction.  From  this  study,  the 
major  contribution  is  the  description  of  an  Eulerian  return  flow  (Hasselman’s  drift)  that 
cancels  the  Stokes  drift,  allowing  for  a  measure  of  the  uncertainty  related  to  the  waves’ 
influence  on  the  currents. 

5.  Yelland  and  Taylor,  1995 

This  is  a  work  on  wind  stress  measurements  from  the  open  ocean.  The  values  of 
the  drag  coefficient  (Co)  for  the  theoretical  approach  made  in  II. D.  are  used  according  to 
this  study’s  conclusions: 

•  For  winds  speeds  at  10  m  (Uio)  from  3  to  6  m/s 

1000Co  =  0.6  +  0.07U10  (5) 

•  For  winds  speeds  at  10  m  from  6  to  26  m/s 

1000Co  =  0.6  +  0.07U10  (6) 

6.  Ardhuin,  Chapron  and  Elfouhaily,  2003 

This  is  a  work  about  the  interaction  between  wind  and  waves.  The  result  of  this 
study  that  most  contributes  to  this  thesis  is  that  for  growing  waves  in  short  fetches,  there 
is  a  wave  induced  stress  that  opposes  the  wind  stress,  consuming  up  to  10%  of  the 
energy.  This  directly  influences  the  transfer  function  and  the  forecast  that,  practically,  is 
only  based  on  actual  currents  and  wind  forecast.  To  remove  this  influence,  a  measure  of 
the  wave  height  difference  due  to  wind  presence  should  be  introduced  -  which  it  self  is  a 
function  of  magnitude  and  direction  of  the  wind  and  of  the  waves  height,  wavelength  and 
direction  and  of  the  fetch. 

7.  Ullman  et  al.,  2003 

These  authors  prepared  a  for  the  United  States  Coast  Guard,  Department  of 
Homeland  Security,  on  the  use  of  CODAR  technology  in  Search  and  Rescue  (SAR) 
planning.  This  work  was  continued  by  the  next  study  referenced. 
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8.  O’Donnall  et  al.,  2005 

This  is  also  a  study  prepared  for  the  United  States  Coast  Guard,  Department  of 
Homeland  Security,  on  the  integration  of  CODAR  and  a  Short-Term  Predictive  System 
(STPS)  surface  current  estimate  into  the  Coast  Guard’s  SAR  Optimal  Planning  System 
(SAROPS).  This  study  is  a  catalyst  and  at  the  same  time  a  measure  of  comparison  for  the 
results  of  this  thesis.  As  described  in  I.B.2.,  its  goals  and  approach  were  the  same,  but 
used  more  limited  data  and  a  different  mathematical  analysis.  One  of  its  conclusions 
mentioned  the  necessity  of  a  work  such  as  this  thesis.  This  thesis  results  ended  up  being 
similar. 


9.  Garfield,  Paduan  and  Ohlmann,  2007 

A  report  prepared  for  the  Coastal  Response  Research  Center,  this  is  a  study 
funded  by  NOAA  and  the  University  of  New  Hampshire.  One  of  the  focuses  of  this 
report  was  to  test  the  accuracy  of  the  Naval  Postgraduate  School  model,  and  it  pointed  to 
the  necessity  of  studying  the  other  interactions,  starting  with  the  wind.  For  this  purpose,  a 
comparison  is  made  between  CODAR  forecasted  positions  and  CODAR  reading 
positions  of  simulated  drifters  after  24h.  Its  conclusions  mention  that  “one  environmental 
factor  that  is  expected  to  play  a  role  is  wind,”  despite  inconclusive  results.  The  study  also 
mentions  that  there  is  “some  connection  between  strong  currents  (...)  and  large  RMS 
separation.” 

This  thesis  is  an  active  part  of  the  continuation  of  this  work. 

10.  Ardhuin  et  al.,  October  2008 

This  study  came  to  the  attention  of  the  author  during  the  writing  stage  of  this 
thesis.  Coincidentally,  it  concerned  exactly  the  focus  of  section  II:  the  inducing  of 
currents  by  the  wind.  Its  conclusions  bolster  the  ones  reached  in  this  thesis,  having  points 
in  common  such  as  the  magnitude  of  the  wind  driven  current  measured  around  2%  and 
the  angle  between  the  velocity  vectors  observed  from  10  to  40  degrees,  based  on  a  two- 
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year  time-series  analysis.  Moreover,  he  study  goes  further  in  theorizing  about  the  wave’s 
influence,  one  of  this  thesis’s  proposals  for  necessary  future  work  ,  as  pointed  to  by 
II.A.4.  and  II.A.6. 

B.  DATA 

As  mentioned  earlier,  this  work  is  one  of  the  first  able  to  study  two  long  reliable 
time  series,  with  continuous  one  hour  measures  for  more  than  one  year.  All  data  analyzed 
was  gathered  from  June  1,  2006  to  December  31,  2007.  Figure  2  shows  the  time  intervals 
for  each  station.  This  allows  a  wide  spectrum  of  possibilities  and  test  scenarios  in  space 
and  time,  as  well  as  providing  outputs  of  the  time  series  analysis  with  a  very  good 
resolution,  maintaining  a  continuous  information  flow  with  a  high  confidence  level. 
Current  data  is  gathered  from  10  CODAR  stations,  with  calculations  of  currents  for  five 
NDBC  buoys  positions,  from  where  the  wind  data  are  obtained  (Figure  3). 

Time  Series  intervals  of  matched  Wind  and  Current  data 

Ml 

mo  y  - 

HMB - 

MRY 1 

I  » _ I _ I  I _ I _ I 

Jul06  Oct06  Jan07  Apr07  Jul07  Oct07  Jan08 

Figure  2.  Matching  data  of  Wind  and  Current  -  Each  jump  in  the  line  corresponds  to 
an  interval  with  no  data. 
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Figure  3.  Current  and  Wind  reading  stations  -  Dots  represent  one  example  hourly 
coverage  of  the  Pescadero  CODAR  station.  (From:  ArcGIS.  Coastline  from 
NOAA,  National  Ocean  Service,  1994.) 

1.  Winds 

Winds  were  obtained  from  NOAA’s  National  Data  Buoy  Center  (NDBC;  Figure 
4,  Table  1),  and  from  the  Monterey  Bay  Aquarium  Research  Institute  (MBARI;  Figure  5, 
Table  2)  moored  buoys. 
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a . 


NDBC  Buoys 


Figure  4.  NDBC  type  buoy.  (From:  National  Data  Buoy  Center, 
http://www.ndbc.noaa.gov,  October  2008). 


Data  from  the  following  NDBC  buoys  /  stations  is  used: 


# 

Name 

Distance  to 
shore 

Type 

Position 

Elevation 

Anemometer 

height 

Water 

depth 

Watch 
circle  radius 

46012 

Half  Moon  Bay 
(HMB) 

4  NM 

3  m 
discus 

37.361N 

122.881W 

sea 

level 

5  m 

213.1  m 

170  yds 

46026 

San  Francisco 
(SFR) 

18  NM 

3  m 
discus 

37.759N 

122.833W 

sea 

level 

5  m 

52.1  m 

127  yds 

46042 

Monterey 

(MRY) 

27  NM 

3  m 
discus 

36.753N 

122.423W 

sea 

level 

5  m 

2115m 

2175  yds 

Table  1.  NDBC  Buoys  (From:  National  Data  Buoy  Center, 

http://www.ndbc.noaa.gov,  October  2008.) 
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b.  MBARI  Buoys 


Figure  5.  MBARI  MO  and  Ml  buoys.  (From:  http://www.mbari.org.  October  2008). 
Data  from  the  following  MBARI  buoys  /  stations  is  gathered  from: 


Name 

Position 

Elevation 

Anemometer 

height 

Water 

depth 

MO 

6.835N 

121.899W 

sea 

level 

4  m 

70  m 

Ml 

36.750N 

122. 020 W 

sea 

level 

4  m 

1000  m 

Table  2.  MBARI  Buoys  (From:  Monterey  Bay  Aquarium  Research  Institute, 
http://www.mbari.org,  October  2008;  Depths  (From:  ICON, 

http://www.oc.nps.edu/~icon/moorings,  November  2008). 

Notice  that  since  MO  and  Ml  are  well  inside  the  Monterey  Bay  (Figure  5), 
their  data  differ  from  the  other’s  positions  because  they  are  under  the  control  of  the 
typical  Bay  circulation  and  winds,  of  high  topographic  and  land  influence. 

c.  Resolution  and  Accuracy 

The  resolution  of  the  anemometers  is  unknown  but  the  readings  are  valued 
to  0.0001  m/s.  The  accuracy  is  also  unknown,  but  considering  the  readings’  detail, 
variance  is  considered  to  be  small. 
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d.  Handling 

In  order  to  enable  treatment  and  comparison  among  different  equipment 
readings,  all  the  data  must  be  formatted  into  a  similar  shape.  Therefore,  all  of  the  wind 
time  series  are: 

•  Hourly  averages  of  the  available  10  minute  readings,  from  the  exact  hour  - 
30  min  to  +  30  min; 

•  Arranged  in  the  form  of  a  three  dimensional  vector  of  time  (in  hours),  u 
velocity  (seasonal,  in  m/s)  and  v  velocity  (meridional,  in  m/s); 

•  Interpolated  for  gaps  under  6  hours; 

•  arranged  so  that  seasonality,  monthly,  synoptic  and  daily  influences  are 
studied  by  subdivisions  of  the  series  in  intervals  of  three  months,  one 
month,  six  days  and  one  day. 

•  calibrated  such  that  spatial  influence  is  studied  by  the  different  locations. 

e.  Other  Series  Created 

A  different  time-series  is  also  created  for  the  wind  readings  of  the 
Monterey  and  Half  Moon  Bay  NDBC  buoys,  where  the  y  axis  is  rotated  for  the  winds’ 
largest  variance  azimuth,  maintaining  the  x  axis  in  perpendicular  alignment.  This  will 
serve  to  interpret  the  angle  between  the  two  major  variances  of  wind  and  current, 
extracting  conclusions  on  interaction  both  for  rotated  axis  for  wind  biggest  variation  only 
and  for  wind  data  rotated  for  the  wind  biggest  variation  and  current  to  the  current  biggest 
variation. 

2.  Currents 

Currents  were  obtained  from  all  the  in-range  CODAR  stations,  for  all  the 
positions  of  the  wind  stations  mentioned  in  the  previous  section  from  the  buoys  positions. 
A  photograph  of  a  typical  CODAR  transmit  and  receive  antenna  configuration  is  shown 
in  Figure  6.  Despite  the  increasing  number  of  CODAR  stations  around  the  world,  there 
have  not  been  many  opportunities  to  operate  for  a  long  period  with  so  many  stations  in  a 
single  study,  increasing  the  number  of  radial  readings  for  each  position.  A  short 
explanation  of  how  CODAR  technology  appeared  and  operates  is  given  in  paragraph  a. 
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fl. 


CODAR  Technology 


Figure  6.  Point  Sur  CODAR  station  -  The  receiver  antenna  is  closer  and  the 
transmitter  further  down.  (From:  CODAR,  http://www.codaros.com, 
October  2008). 

CODAR  consists  of  radar  in  the  HF  band,  which  is  the  band  that 
corresponds  to  wavelengths  of  the  same  magnitude  of  the  sea  waves.  The  first 
correlations  between  radar  and  waves  appeared  accidentally  at  the  end  of  the  Second 
World  War,  when  echoes  appeared  in  near  shore  radar  stations.  In  1955,  “Crombie 
discovers  experimentally  that  Bragg  scatter  produces  strong  FIF  sea  echo”  (CODAR, 
http://www.codaros.com,  October  2008). 


CODAR  uses  the  Bragg  Scattering  phenomenon  (Figures  7  and  8)  to 
measure  the  wavelength  of  the  sea  waves  that  are  approaching  or  going  away,  in  a  radial 
direction,  to  each  station.  Having  the  wavelength,  we  calculate  the  speed  a  wave  should 
have  from  the  deep  water  dispersion  relation.  The  difference  of  where  the  received 
frequency  peak  is  and  where  it  should  be  in  that  particular  wavelength  is  a  measure  of  the 
radial  surface  current  by  (Doppler,  1803-1853): 


/'  =  /„ 


1  + 


(7) 


19 


/  ’  being  the  new  peak  frequency,  /0  the  expected  peak  with  no  current,  v 
the  radial  current  velocity  and  c0  the  theoretical  velocity  of  a  wave  with  that  wavelength 

in  deep  water.  With  two  of  these  measurements,  if  properly  intercepting,  we  obtain  the 
two  components  of  the  velocity  vector. 


Figure  7.  Bragg  scattering  phenomenon  -  The  biggest  return  in  scattered  energy  is  Vi 
X  (wavelength)  of  the  sea  wave.  (From:  Martin,  2004  in  “An  introduction  to 
Ocean  Remote  Sensing,”  Cambridge.) 


20 


Figure  8.  Bragg  scattering  phenomenon.  Received  power  by  frequency  -  The 
difference  of  where  the  peak  is  and  where  it  should  be  for  that  particular  X  , 
is  proportional  to  the  radial  current  velocity.  (From:  Paduan  &  Graber, 
1997.) 

For  the  frequencies  used  in  these  stations,  the  measurements  are  related  to 
the  first  one  meter  of  the  water  column;  horizontally,  the  result  for  a  single  observing  bin 
is  based  on  backscatter  from,  approximately,  3x3  km  of  surface  area. 

b.  CODAR  Stations 

All  stations  but  the  one  off  Moss  Landing  consisted  of  a  standard  1 1.5-14 
MFIz  apparatus;  Moss  Landing  operated  in  the  range  24-27  MHz.  Their  approximate 
positions  (in  Geographic  Coordinate  System)  are  given  in  Table  3. 
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CODAR  station 

Position 

Commonwheel 

37.91 1833N 

122.728167W 

Fort  Funs  ton 

37.712050N 

122.5013W 

Montara 

37.533717N 

122.519217W 

Pescadero 

37.252483N 

122.416567W 

Santa  Cruz 

36.9492 17N 

122.0661W 

Moss  Landing 

36.803667N 

121.787917W 

Point  Pinus 

36.636783N 

121.9536W 

Naval  Postgraduate  School 

36.603333N 

121.87195W 

Granite  Canyon 

36.439450N 

121.922217W 

Point  Sur 

36.305767N 

121.90125W 

Table  3.  CODAR  stations  (From:  http://www.cencalcurrents.com,  September  2008.) 


22 


c.  Resolution  and  Accuracy 

The  resolution  for  these  stations  is  0.04  m/s  (from 
http://www.codaros.com.  accessed  in  October  2008).  The  most  recent  studies  show  accuracy 
from  8  to  18  cm/s  (from  Wright,  2008,  NPS  Masters  Thesis). 

d.  Handling 

In  order  to  be  able  to  work  with  these  time  series  all  the  data  have  to  be 
arranged  in  a  similar  fashion.  All  of  the  current  time  series  are: 

•  Hourly  averages  from  the  exact  hour  -30  min  to  +  30  min,  interpolated 
from  the  nearest  valued  positions  to  the  buoys  position.  This  is  the  way  the 
data  is  presented  by  the  CODAR  system. 

•  Arranged  in  the  form  of  a  three  dimensional  vector  of  time  (in  hours),  u 
velocity  (seasonal,  in  m/s)  and  v  velocity  (meridional,  in  m/s).  This 
arrangement  allows  easy  handling. 

•  Interpolated  for  gaps  under  6  hours.  This  allows  a  good  continuity  and 
guarantees  no  hiding  of  important  features  in  the  changing  of  the  current. 

e.  Other  Series  Created 

Other  time-series  were  also  created  for  the  current  readings  of  the  position 
of  the  Monterey  and  Half  Moon  Bay  NDBC  buoys: 

•  Where  the  v  axis  is  rotated  for  the  wind’s  largest  variance  azimuth, 
maintaining  a  perpendicular  x  axis; 

•  Where  the  y  axis  is  rotated  for  the  current’s  largest  variance  azimuth, 
maintaining  a  perpendicular  x  axis.  (Comparison  of  the  results  from  steps 
(1)  and  (2)  also  provides  information  about  the  most  common  angle 
between  surface  current  and  wind.) 

•  All  the  series  created  for  the  Monterey  buoy  have  two  more  associated 
scenarios,  where  the  tide  influence  from  the  current  is  tentatively  removed. 
For  this  purpose  the  software  t  tide  is  used  (Pawlowicz,  Beardsley  and 
Lentz,  2002  with  a  Cook,  2008,  alteration  to  the  software).  In  the  first,  all 
of  the  tide  constituents  available  in  the  software  are  used,  once  this  time- 
series  is  long  enough  (in  this  case,  one  and  a  half  years).  In  the  second 
one,  all  but  the  SI  constituents  are  used,  in  order  to  avoid  removing  the 
sea  breeze  influence.  In  both  we  have  to  take  into  account  that  the  energy 
removed  may  not  be  just  tide  related,  as  occurs  in  periods  of  the  order  of 
the  rotation  of  the  planet’s  rotation  frequency. 
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•  Seasonality,  monthly,  synoptic  and  daily  influences  are  studied  by 
subdivisions  of  the  series  in  intervals  of  three  months,  one  month,  six  days 
and  one  day. 

•  Spatial  influence  is  studied  by  the  different  locations. 

3.  Both  Wind  and  Current 

The  time  series  of  both  wind  and  current  are  then  cut  in  order  to  create  time 
intervals  where  they  both  exist,  to  make  comparison  possible.  The  resulting  matched 
intervals  are  displayed  in  Figure  2. 

From  now  on  in  this  thesis,  for  reference  simplicity,  when  a  station  is  mentioned 
by  place  name  it  refers  to  the  buoy  position;  e.g.,  “Monterey  Station”  means  the  data 
obtained  for  the  position  of  the  46042  NDBC  buoy.  Because  the  Monterey  station  has  a 
broader  continuous  time  series,  it  allows  us  to  extrapolate  more  results  with  a  better 
resolution.  That  particular  location  is  where  most  studies  and  essays  are  made,  the  others 
serving  as  comparisons. 

C.  EMPIRICAL  APPROACH 

Two  independent  paths  are  taken.  The  first  starts  with  real  data  and  tries  to 
extrapolate  relations  from  it.  This  approach  relies  on  pure  Time  Series  Analysis.  The 
second,  described  in  II  D.,  is  based  on  theoretical  knowledge  and  serves  as  comparison  to 
the  empirical  results. 

All  of  the  following  methods  are  applied  in  all  of  the  time  series  described  in  B. 

1.  Principal  Axis  (Biggest  Variance) 

The  objective  is  to  have  a  measure  of  the  axis  which  experiences  the  largest 
variance  and  its  angle  of  measurement,  for  both  wind  and  current.  The  result  is  plotted  as 
Figure  9  and  as  a  table  in  Table  4,  where  the  two  major  semi-axes  in  the  along  vector 
direction  are  shown  in  full  (using  Bahr,  2007,  software): 
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Elipses  axis  from  01 -Jul-2006  +  1.5  years  :  alfa=-58.5547deg 


U 

Figure  9.  Principal  axis  analysis  -  Values  correspond  to  variance  and  alfa  is  the  angle 
between  the  two  full  lines.  (Monterey  station.) 

The  table  shows  values  for  each  term: 

Principal  axis  for  W/C:  from  01-Jul-2006  +  1.5  years 


theta 

Maj. 

min. 

wind 

-30.4301 

4.8623 

1.9311 

current 

28.1246 

0.13539 

0.13232 

Table  4.  Principal  axis  analysis.  -  Values  correspond  to  magnitudes  and  angle 
between  the  major  axis  and  north.  (Monterey  station.) 


The  results  of  all  the  time  series  show  that  magnitudes  of  the  variance  keep  within 
the  same  order.  The  angles  vary  more,  getting  dispersed  as  we  reduce  the  time  series 
length.  The  average  angle  between  the  two  major  variance  directions  is  around  50 
degrees  in  the  longer  time-series  as  shown  in  Table  5. 
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Examples  of  results  from  the  different  stations  in  different  periods 


Station 

Period 

Angle  between  major  axis  (°) 

MRY 

01JUF06  ..  31DEC07 

-59 

HMB 

01JUF06  ..  23NOV06 

-14 

20FEB07  ..  12JUE07 

-8 

04SEP07  ..  31DEC07 

-10 

SFR 

29NOV06  ..  07JAN07 

-2 

01MAR07  ..  26MAY07 

-15 

06SEP07  ..  31DEC07 

-19 

MO 

1 1OCT06  ..  28MAR07 

-39 

29MAR07  ..  25JUN07 

+21 

25JUN07  ..  31DEC07 

+25 

Ml 

16NOV06  ..  05JAN07 

-43 

24APR07  ..  24JUN07 

+  1 

25JUN07  ..  30AUG07 

-6 

Table  5.  Principal  axis  analysis.  -  Space  and  time  variability  (wind  -  current.) 

Despite  the  tendency  for  the  longer  time  series,  the  high  variability  does  not  allow 
a  good  estimation  of  a  typical  angle  between  wind  and  current,  demonstrating  that  it  may 
be  more  a  function  of  time  and  place. 

2.  Complex  Correlation  (Kundu,  1975) 

Complex  correlation  gives  a  measure  of  correlation  between  two  complex  time 
series.  The  magnitude  of  the  result  shifts  between  0  and  1,  1  being  if  the  signal  is  the 
same.  The  average  phase  angle  gives  an  idea  of  the  physical  angle  between  the  two 
phenomena.  Results  for  the  Monterey  station  are  plotted  on  Figure  10,  and  results  for  all 
stations  are  tabulated  in  Table  6. 
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Complex  Correlation  W/C  MRY  JUL06  to  DEC07 


0.8 
0.6 
0.4 
0.2 
jS  0 
-0.2 
-0.4 
-0.6 
-0.8 
-1 

-1  -0.8  -0.6  -0.4  -0.2  0  0.2  0.4  0.6  0.8  1 

Real 

Figure  10.  Complex  correlation  -  Modulus  is  a  measure  of  correlation  and  alpha  is  the 
physical  angle.  (Monterey  station.) 

The  magnitudes  vary  from  0.2  to  0.7.  The  angles  tend  to  remain  on  the  second 
quadrant  around  35  degrees,  gaining  variance  with  the  reducing  of  the  length  of  the  time 
series.  These  results  mean  that  they  are  correlated,  and  there  is  in  average  35  degrees 
between  the  directions  of  the  two.  This  is  expected  because  the  angle  is  supposed  to 
increase  in  time,  until  Ekman’s  45  degrees.  The  results  for  M0  buoy  present  the  highest 
variability  but  also  the  smaller  correlation,  probably  due  to  the  variability  of  the  winds 
and  currents  inside  the  Bay.  An  attempt  was  made  to  verify  if  the  correlation  would 
increase  when  lagging  the  Current  time  series,  but  it  decreased,  pointing  to  zero  lag. 


Examples  of  results  from  the  different  stations  in  different  periods 


Station 

Period 

Magnitude 

Angle  (°) 

MRY 

01JUL06  ..  31DEC07 

0.49 

-32 

HMB 

01JUL06  ..  23NOV06 

0.31 

-56 

20FEB07  ..  12JUL07 

0.66 

-47 

04SEP07  ..  31DEC07 

0.48 

-56 

SFR 

29NOV06  ..  07JAN07 

0.66 

-1 

01MAR07  ..  26MAY07 

0.81 

-16 

06SEP07  ..  31DEC07 

0.62 

-23 

i - 1 -  ,  i  i 

- mod=0.48542  alfa=— 31 .7407 
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Station 

Period 

Magnitude 

Angle  (°) 

MO 

1 1OCT06  ..  28MAR07 

0.34 

-26 

29MAR07  ..  25JUN07 

0.27 

-83 

25JUN07  ..  31DEC07 

0.17 

-45 

Ml 

16NOV06  ..  05JAN07 

0.5 

-34 

24APR07  ..  24JUN07 

0.67 

-17 

25JUN07  ..  30AUG07 

0.69 

-13 

Table  6.  Complex  correlation.  -  Space  and  time  variability. 

3.  Scatter  Comparison 

Scatter  plots  give  an  idea  of  correlation  and  pennit  a  mathematical  best-fit 
expression  (e.g.,  Figure  11,  Figure  12,  Figure  13).  A  perfect  correlation  happens  when  all 
the  pairs  of  values  (points)  of  the  two  time  series  design  a  mathematical  expression.  This 
will  allow  elaboration  of  the  best  transfer  function  relation  in  terms  of  magnitude.  The 
one  presented  here  considers  the  full  energy  spectrum,  and  the  best  fit  is  accomplished 
with  a  quadratic  expression. 

This  analysis  can  also  be  done  by  frequency  bands  if  the  choice  is  to  transfer  by 
frequency.  As  we  will  see  later,  the  degree  of  added  complexity  for  such  an  approach  will 
not  increase  significantly  the  results  of  the  first  one. 
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|Wind|  (m/s) 


Figure  13.  Scattering  diagram  for  Modulus  and  best  fit  expression.  (Monterey  station, 
Jul  06  to  Dec  07.) 

After  applying  it  to  all  the  series,  including  the  rotated  ones,  the  best  proxy  for  the 
typical  expression  relating  near  surface  wind  to  surface  current,  in  terms  of  magnitude  of 
the  vector,  is: 

C  =  0.0  ISU5  +  0.03  (8) 

where  C  is  the  surface  current  in  m/s  and  U5  is  the  wind  at  5  m  height,  also  in  m/s. 

The  quadratic  term  shown  on  the  figures  is  always  positive  and  negligible,  with 
some  exceptions  in  the  modulus  when  it  goes  negative.  The  results  presented  here  are 
restricted  to  the  linear  tenns  in  the  fit.  Those  values  multiplying  by  U5  vary  plus  or  minus 
0.02.  But,  overall,  the  results  show  surface  current  magnitudes  are  very  close  to  2%  of  the 
wind  speed  magnitudes  (Table  7).  The  choice  of  0.018  is  supported  by  the  error  analysis 
and  tuning  done  in  the  forecast  section.  The  last  constant  varies  plus  or  minus  0.05.  This 
variation  tends  to  be  smaller  in  winter  and  fall  and  in  the  fits  based  on  longer  time  series. 
This  expression  will  be  tuned  in  the  predicting  algorithm  in  Section  III. 
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Examples  of  results  from  the  different  stations  in  different  periods 


Station 

Period 

U 

V 

Modulus 

MRY 

01JUL06  ..  31DEC07 

0.006X-0.028 

0.018X+0.018 

0.007X+0.128 

HMB 

01JUL06  ..  23NOV06 

0.001X-0.063 

0.02 1X+0. 096 

-0.016X+0.217 

20FEB07  ..  12JUL07 

0.000X-0.047 

0.017X-0.068 

0.001X+0.178 

04SEP07  ..  31DEC07 

-0.006X-0.045 

0.018X-0.013 

0.007X+0.167 

SFR 

29NOV06  ..  07JAN07 

0.013X-0.031 

0.022X+0.037 

-0.016X+0.162 

01MAR07  ..  26MAY07 

0.013X+0.000 

0.018X-0.033 

0.006X+0.175 

06SEP07  ..  31DEC07 

0.007X-0.012 

0.023X+0.004 

0.013X+0.106 

MO 

11OCT06  ..  28MAR07 

0.005X-0.014 

0.007X+0.001 

0.001X+0.099 

29MAR07  ..  25JUN07 

0.001X-0.089 

0.004X+0.032 

0.006X+0.1 14 

25JUN07  ..  31DEC07 

0.004X-0.039 

0.006X+0.017 

0.006X+0.104 

Ml 

16NOV06  ..  05JAN07 

0.008X-0.039 

0.016X+0.041 

0.007X+0.1 1 1 

24APR07  ..  24JUN07 

-0.009X+0.025 

0.008X-0.028 

-0.005X+0.149 

25JUN07  ..  30AUG07 

0.005X+0.001 

0.010X+0.007 

-0.015X+0.177 

Table  7.  Best  fit  expression.  -  Space  and  time  variability. 

4.  Auto  Correlation  (McMahan  and  Wyland,  NPS  Lab  Software,  2007) 

Autocorrelation  of  a  given  time  series  allows  for  a  better  understanding  of  the 


time  scales  involved.  Each  oscillation  represents  the  time  scale  of  an  active  agent. 
Examples  are  shown  for  wind  and  current  in  Figure  14  and  Figure  15,  respectively. 
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auto  correlation  Wind  U  from  01Uul-2006  +  24h 
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auto  correlation  Wind  V  from  01  -Jul-2006  +  24h 
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auto  correlation  Wind  U  from  01Uul-2006  +  24h 


auto  correlation  Wind  V  from  OIUul-2006  +  24h 


Figure  14.  Autocorrelation  for  Wind  -  The  bottom  is  a  zoom  of  the  first  positive  2000h 
~  3  Months.  Daily,  synoptic  and  monthly  oscillations  are  present.  (Monterey 
station,  Jul  06  to  Dec  07.) 
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Concerning  the  wind,  the  influence  of  the  sea  breeze  is  easily  seen  on  the  U 
component  in  the  small  oscillations  of  24  h.  On  both,  we  can  also  see  an  oscillation 
around  the  6  days  (synoptic  scale).  In  U,  monthly  and  a  seasonal  oscillations  are  also 
present. 
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auto  correlation  Current  V  from  01Uul-2006  +  24h 
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auto  correlation  Current  U  from  01  -Jul-2006  +  24h 


auto  correlation  Current  V  from  01  ^Jul-2006  +  24h 


Figure  15.  Autocorrelation  for  Current  -  The  bottom  is  a  zoom  of  the  first  positive 
2000h  ~  3  Months.  Daily,  synoptic  and  monthly  oscillations  are  present. 
(Monterey  station,  Jul  06  to  Dec  07.) 
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For  the  current,  daily,  synoptic  and  monthly  oscillations  are  noticed  as  well.  The 
tidal  influence  can  be  easily  seen  from  the  difference  between  Figure  15  (simple  current 
measurement)  and  Figure  16  (after  tide  removal).  Daily,  synoptic  and  monthly 
oscillations  are  present,  but  the  daily  one  is  much  less  energetic.  All  the  daily  oscillations 
tend  to  be  bigger  in  the  summer.  There  is  no  important  change  looking  between  stations, 
except  in  the  magnitude  of  the  oscillations. 


Tides  out  -  auto  correlation  Current  V  from  01-Jul-2006  +  24h 


Tides  out  -  auto  correlation  Current  U  from  01^Jul-2006  +  24h 


Tides  out  -  auto  correlation  Current  V  from  01-Jul-2006  +  24h 
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Figure  16.  Autocorrelation  for  Current  with  tides  removed  -  The  bottom  is  a  zoom  of 
the  first  positive  2000h  ~  3  Months.  Daily,  synoptic  and  monthly 
oscillations  are  present.  (Monterey  station,  Jul  06  to  Dec  07.) 
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5. 


Power  Spectra  and  Straight  Coherence 


The  power  spectra  of  wind  and  current  give  the  distribution  of  energy  versus 
frequency  (Figures  17  and  18).  Based  on  this,  coherence  gives  the  significance  of  the 
relation  of  the  signals  for  each  frequency  and  the  phase  difference  between  them  (Figure 


Figure  17.  Power  spectra  of  Wind  and  Current  -  Dashed  at  95%  confidence.  (Monterey 
station,  Jul  06  to  Dec  07.) 

The  peaks  for  the  daily  and  12h  periods  are  easily  spotted  in  both  wind  and 
current.  Winds  are  especially  energetic  in  the  sea  breeze  period  while  the  currents  show 
almost  equivalence  on  both,  maybe  because  of  the  previously  mentioned  wind  and  tide 
correlations.  After  the  tide  signal  removal  attempt,  both  current  peaks  are  slightly 
smaller,  representing  less  energy.  When  removing  tides  by  best  fit  of  frequencies,  we  are 
removing  also  part  of  the  wind  influence  which  shares  those  same  frequencies  -  the  sea 
breeze. 
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Tides  out  -  Gxmean  Gymean  MRY 


Figure  18.  Power  spectra  of  Wind  and  Current  with  tides  removed  -  Dashed  at  95% 
confidence  (Monterey  station,  Jul  06  to  Dec  07.) 

Adding  the  variances  both  with  and  without  tides,  we  notice  the  slight  decrease  in 
energy  that  is  removed.  The  results  are  summarized  in  Table  8. 


Difference  of  with  and  without  tides 


wind  variance 

19.5694 

tides  in  -  current  variance 

0.016719 

tides  out  -  current  variance 

0.015593 

Table  8.  Added  variances  ~  energy.  (Monterey  station.) 

These  series  are  significantly  coherent  for  periods  longer  than  a  day.  During  this 
interval,  the  phase  augments  from  near  zero  to  around  60°.  In  the  shorter  time  series  the 
coherence  is  consistent  only  after  40  h,  and  the  phase  tends  not  to  exceed  50°.  When  done 
by  components  (u  and  v)  it  is  coherent  after  approximately  12  h. 
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Figure  19.  Coherence  and  Phase  between  wind  and  current  -  Dashed  at  95% 
confidence.  (Monterey  station,  Jul  06  to  Dec  07.) 


6.  Rotary  Spectra  (Gonnela,  1972  and  Mooers,  1972,  Bahr,  2007 
Software) 


Rotary  spectra  subdivides  a  signal  into  two  rotating  components,  one  clockwise 
and  the  other  counter-clockwise.  This  allows  one  to  analyze  the  coupling  of  two  signals 
in  terms  of  their  rotating  behavior,  giving  independence  from  the  axis  choice.  Results  for 
wind  and  current  for  the  Monterey  station  are  shown  in  Figure  20  and  Figure  21, 
respectively. 


Figure  20.  Rotary  components  for  wind  -  Notice  the  two  peaks  for  12  and  24h. 
(Monterey  station,  Jul  06  to  Dec  07.) 
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Figure  21.  Rotary  components  for  current,  tides  in  and  out  -  Notice  a  small  decrease  in 
energy.  (Monterey  station,  Jul  06  to  Dec  07.) 

The  two  diurnal  and  semi-diurnal  peaks  are  distinct  in  both  wind  and  current, 
reflecting  the  importance  of  such  frequencies.  Notice  that  the  clockwise  component  is 
bigger,  the  peaks  being  more  similar  to  the  24h  point  in  the  wind  and  the  12h  in  the 
current.  The  current  peaks  are  only  a  little  smaller  after  the  tides’  removal. 

The  “rotary  coefficient”  is  a  measure  of  the  purity  of  the  signal  in  terms  of 
rotation:  a  fully  rotating  signal  is  1,  a  non-rotating  is  0.  The  rotary  coefficient  for  the 
wind  and  current  records  at  the  Monterey  station  are  shown  in  Figure  22  and  Figure  23, 
respectively.  Distinctively,  for  periods  under  12h,  the  current  is  more  rotational  than  the 
wind,  which  only  shows  two  small  peaks  in  the  frequencies  of  interest. 
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Tides  in  -  Rotary  Coefficient  for  W  MRY.  From  01-Jul-2006  +  1.5  years 


Figure  22.  Rotary  coefficient  for  Wind  -  A  fully  rotational  movement  would  be  valued 
1.  (Monterey  station,  Jul  06  to  Dec  07.) 


Tides  in  -  Rotary  Coefficient  for  C  MRY.  From  01-Jul-2006  +  1.5  years 


Tides  out  -  Rotary  Coefficient  for  C  MRY.  From  01-Jul-2006  +  1.5  years 


Figure  23.  Rotary  coefficient  for  Current,  tides  in  and  out.  (Monterey  station,  Jul  06  to 
Dec  07.) 
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7.  Vector  Cross  Spectral  Analysis  (Mooers,  1972  and  McNally  et  al, 
1989,  After  Jessen  and  Bahr,  2007  software) 

Vector  cross  spectral  analysis  is  a  tool  to  check  the  coherence  between  two  rotary 
signals.  Plugging  in  wind  and  current,  the  coherence  results  for  the  Monterey  station  are 
shown  in  Figure  24. 


Phase 
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Tides  in  -  Coherence  for  MRY  from  01-Jul-2006  +  1.5  years 


Figure  24.  Complex  Coherence  and  Phase  -  Lower  image  is  a  zoom  of  the  coherent 
part.  (Monterey  station,  Jul  06  to  Dec  07.) 

The  application  of  this  method  allows  one  to  calculate,  besides  the  coherence,  the 
angular  shift  in  time.  Table  9  shows  the  average  results  for  the  coherent-positive  part. 


40 


Complex  Coherence  results 


Time  (h) 

Angle  O 

20 

30-35 

30 

40-45 

60 

45-60 

Table  9.  Wind  and  Current  shifting  of  angles  in  time.  (Averages  for  all  the  stations 
and  periods.) 


After  17h  almost  all  the  series  are  coherent  and  have  a  much  slower  variation.  For 
periods  between  0  and  17h  the  results  were  inconclusive.  This  treatment  is  done  also  to 
the  velocity  square,  to  search  for  a  difference  between  pure  wind  velocity  and  wind  stress 
(function  of  the  velocity  square).  No  significant  difference  was  found  between  the 
coherence  results  for  wind  or  wind  stress.  For  the  de-tided  series,  the  relation  keeps  its 
results  except  during  the  24  and  12h  periods,  where  very  narrow  positive  and  negative 
peaks  appear  (Figure  25),  probably  a  result  of  energy  removal  in  those  frequencies. 


Tides  out  -  Coherence  for  MRY  from  01-Jul~2006  +  1.5  years 


_0.4  -0.3  -0.2  -0.1  0  0.1  0.2  0.3  0.4 

f  (CPH) 


41 


Tides  out  but  SI  -  Coherence  for  MRY  from  01-Jul-2006  +  1.5  years 
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Figure  25.  Complex  coherence  and  Phase,  tide  removed  -  Top,  all  constituents. 

Bottom,  all  constituents  removed  but  SI.  (Monterey  station,  Jul  06  to  Dec 
07.) 

Comparing  results  for  other  stations  (not  shown),  the  M0  buoy  is  only  coherent 
after  24h  at  a  50°  phase.  The  same  can  be  said  for  Ml  in  the  summer. 

The  empirical  approach  gave  us  estimation  of  the  interaction  between  the  wind 
and  the  surface  current  in  its  magnitudes  and  phase  shifts.  This  allows  building  a  transfer 
function  that  will  later  be  used  to  forecast  the  current,  having  the  forecasted  winds  in 
consideration  of  various  factors. 

D.  THEORETICAL  APPROACH 

In  this  approach,  theory  is  used  to  compare  with  the  empirical  results. 
Calculations  are  made  using  Ekman’s  works  theory,  and  a  function  is  created 
considering: 


T)x  _  PS-rfJ io  • 

And  at  the  same  time 


(9) 
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.  du 

U  ~  Pa^V  • 


(10) 


So,  adding  the  two  and  rearranging, 

du  =  ^^-dz-Ul02.  (11) 

Po  ^  v 

The  magnitude  of  Uio  is  taken  to  be  approximate  of  U5,  for  it  is  the  top  part  of  the 
exponential  wind  growth  in  height;  the  Yelland  and  Taylor,  1995  drag  coefficient  is  used 
as  function  of  the  wind  speed. 

E.  COMPARISON 

Having  the  empirical  and  the  theoretical  result  transfer  functions,  and  using  Uio  = 
10  m/s  we  obtain: 

1 .  From  the  results  of  the  study  in  equation  (8), 

du=0.21  m/s, 

2.  From  the  theoretical  result  of  ( 1 1), 

a.  Ekman  depth  =  50m  and  Ay=0.01  m  /s  du=0.68  m/s; 

b.  Ekman  depth  =  10m  and  Av=0.01  m2/s  du=0.27  m/s; 

c.  Ekman  depth  =  50m  and  Av=0. 3  m2/s  du=0.02  m/s; 

d.  Ekman  depth  =  10m  and  Ay=0.3  m  /s  du=0.01  m/s. 

This  means  that  the  theoretical  result  is  highly  dependent  on  a  good  estimate  of 
eddy  viscosity,  and  that  using  a  typical  value  (Av~0.01),  Ekman  depth  must  be  much 
shallower.  Using  typical  values,  Ekman  theory  overestimates  the  change  in  the  surface 
current  velocity. 

F.  CONCLUSIONS 

Having  studied  the  wind  /  surface  current  interaction  in  different  perspectives, 
times  and  spaces,  the  relationship  found  can  be  described  as  follows. 
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1. 


Transfer  Function 


a.  Magnitude 

The  surface  current  generated  is  around  2%  of  the  wind  speed. 

b.  Angle 

The  angle  measurement  increases  in  an  unknown  way  until  around  17h. 
From  then  on,  the  angles  vary  from  35°  to  60°,  according  to  Table  9. 

2.  Notes 

Despite  these  results,  the  scatter  plots  show  a  large  dispersion,  which  means  that 
there  is  a  high  variability. 

The  unknown  growth  until  17h  includes  wind  wave  interaction.  It  includes  also 
components  of  transfer  from  0  to  35°  of  some  magnitude. 

In  this  study,  rotating  the  axis  to  biggest  variance  or  long-shore,  cross-shore,  does 
not  bring  any  significant  change,  probably  because  of  the  meridional  characteristic  of  the 
Californian  coast.  Nevertheless,  when  using  both  wind  and  current  axis  to  their  respective 
biggest  variances,  the  component  minor  wind  axis  versus  minor  current  axis  does  tend  to 
zero. 

The  tidal  removal  does  not  bring  major  changes  to  the  transfer  function  as  a  bulk. 
If  transferred  by  frequency,  care  has  to  be  taken  especially  in  the  inertial  periods  and  its 
subdivisions. 
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III.  THE  FORECAST 


A.  FORECAST  WITH  CODAR 

Forecasting  with  a  remote  sensing  means  like  CODAR  has  the  enonnous 
advantages  of  independence  and  timing— independence  because  it  does  not  rely  on  any 
other  system  to  operate,  and  timing  because  it  is  updated  every  hour.  This  is  more  than 
enough  to  get  good  estimates  for  the  surface  currents  variability,  making  data  available 
when  needed.  Nevertheless,  relying  only  on  previous  surface  currents  readings  from  our 
own  system  has  its  dangers  and  limitations.  As  will  be  shown,  coupling  our  readings  with 
those  of  other  systems  yields  more  reliable  and  accurate  results. 

1.  Reliability  of  a  Forecast 

CODAR  has  proved  to  be  a  very  reliable  system,  so  far.  Nevertheless,  its  results 
may  be  misleading.  Once  it  uses  electromagnetic  radiation  to  operate,  it  is  subject  to  all 
its  difficulties  and  laws,  as  described  by  the  following  examples. 

•  If  an  object  with  electromagnetic  potential  and  length  comparable  to  the 
transmitted  wavelength  (approximately  10m)  goes  near  the  antennas,  the 
reception  pattern  changes,  leading  to  deviation  in  the  result; 

•  The  presence  of  a  moving  target,  such  as  a  ship,  may  introduce  an  error 
due  to  strong  backscatter; 

•  A  change  in  the  background  radio  frequency  properties  will  influence  the 
signal  also. 

These  and  other  problems  have  been  addressed  several  times,  and  are  adequately 
explored  in  several  works  (some  of  which  can  be  found  on  http://www.codaros.com).  The 
use  of  a  second  system  allows  us  to  detect  possible  malfunctions  by  comparison,  and 
gives  a  different  set  of  solutions  as  backup.  This  is  always  preferable. 

2.  Accuracy  of  a  Forecast 

The  accuracy  of  the  system  is  dependent  on  the  errors  of  CODAR  (8  to  18  cm/s, 
from  Wright,  2008),  plus  the  ones  associated  with  the  estimation  of  the  next  24h  behavior 
(the  assumptions  and  the  errors  of  the  systems  we  are  building  our  prognostic  on),  plus 
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the  error  of  what  we  cannot  or  do  not  predict.  If  we  introduce  to  our  algorithm  a  good 
description  of  the  future  conditions  and  have  good  relation  functions,  we  will  obtain  a 
more  realistic  result. 

So,  coupling  with  other  systems,  as  upload  or  download,  is  advantageous.  We  can 
obtain  data  from  tides,  winds,  drifters,  current  meters,  and  models— and  can  assimilate  our 
results  back  into  the  models. 

B.  PREVIOUS  ATTEMPTS 

The  previous  attempt  model  is  described  in  I  B.  The  results  have  been  generally 
accurate  for  persistent  weather  conditions,  but  tend  to  get  worse  with  changing  weather 
conditions. 

1.  Problems 

The  model  was  based  only  in  tidal  and  persistence  analysis.  Once  any  of  the 
conditions— wind,  waves,  or  any  other  parameter— changes,  the  persistence  does  not 
follow.  In  weather,  significant  changes  can  happen  in  any  24h  period. 

2.  Results 

As  pointed  in  I.B.,  there  is  still  lack  of  confidence  in  this  model  any  time  there  is  a 
significant  change  in  the  weather  (namely  in  current  velocity  and  possibly  wind).  A 
summary  of  results  for  the  previous  tidal  model  applied  to  the  region  offshore  San 
Francisco  is  shown  in  Figure  26.  Here  a  pixel  grid  of  pseudo  drifters  release  points  was 
used  to  create  predicted  trajectories  for  24h,  which  was  repeated  each  hour  for  the  period 
August  2006  through  July  2007.  The  plot  shows  the  nns  spatial  difference  between  the 
end  points  of  the  predicted  trajectories  compared  with  the  observed  trajectories  one  day 
after.  The  negative  correlation  with  number  of  drifters  (lower  panel)  reflects  the  fact  that 
higher  speed  conditions  result  in  more  pseudo  drifters,  leaving  the  sampling  grid  within  a 
24h  period.  The  direct  effect  of  wind  speed  on  the  spatial  accuracy  of  the  predicted 
trajectories  is  less  clear.  Nonetheless,  wind  (specifically  wind  changes)  is  hypothesized  to 
explain  some  error  in  forecast  currents. 
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Figure  26.  Results  for  RMS  compared  with  wind  and  number  of  pseudo  drifters  -  RMS 
for  24h.  The  number  of  pseudo  drifters  present  may  be  used  as  a  proxy  for 
current  speed.  (Image  from  Garfield  et  al.,  2007.) 

C.  THE  INTRODUCTION  OF  WIND 

One  attempt  to  reduce  forecast  error  is  the  introduction  of  wind  influence.  A  study 
of  the  wind  surface  /  current  interaction  was  made  for  this  thesis  and  is  described  in 
Section  II.  The  objective  was  to  build  a  transfer  function  that  allows  us  to  plug  in  the 
wind  driven  current  into  the  forecast.  Several  approaches  were  considered:  transfer  by 
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frequency,  transfer  by  different  division  of  spatial  and  time  components,  transfer  as  a 
vector,  transfer  with  tide  analysis,  and  transfer  based  on  persistence.  The  approaches  that 
seemed  to  be  more  efficient  for  practical  purposes  (accuracy  versus  simplicity)  narrowed 
to  only  one  wind/current  transfer  function  and  transfer  as  a  vector.  Transfer  by  frequency 
was  not  chosen  because  of  the  relatively  small  magnitude  of  the  end  results.  Tide  removal 
and  simple  persistence  conditions  would  have  yielded  similar  results.  The  models  tested 
were  based  on  the  diagram  in  Figure  27,  regardless  whether  or  not  tides  were  removed. 

The  winds  had  different  methodologies  applied  also  in  tenns  of  what  time  slice 
constituted  future  for  computing  wind  differences  and  whether  or  not  to  use  wind 
difference  compared  with  removing  and  then  reintroducing  the  full  wind-driven  current. 
The  results  were  practically  the  same.  In  this  hindcast  study,  forecast  wind  is  taken  from 
the  actual  readings,  simulating  a  zero  error  wind  forecast.  The  methodology  tested  here 
can,  however,  be  immediately  transferred  to  winds  forecasted  by  an  atmospheric  model 
such  as  the  Navy's  COAMPS  model. 


Figure  27.  Basic  flow  chart  diagram  of  the  forecast  procedure  -  Several  adaptations 
were  attempted. 

1.  With  Tides  Treatment 

This  attempt  corresponds  exactly  to  the  previous  diagram.  The  best  results  were 
obtained  using  the  following  procedure. 

We  are  at  t=0  and  we  want  the  next  24h  surface  currents  forecast,  so  we 
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•  Take  the  last  24h  current  readings; 

•  Remove  the  tides  by  “t  tide”  using  all  the  previous  month’s  currents  to  the 
best  fit.  Conducting  the  detiding  with  using  one  month  of  data  or  just  three 
days  of  data  both  returned  good  results  as  shown  in  Table  10. 

Detidal  associated  errors 


RMS  (m/s) 

3  Days 

1  Month 

U 

0.1087 

0.1073 

V 

0.1027 

0.1145 

modulus 

0.1515 

0.1593 

Table  10.  Detidal  errors  with  3  days  and  1  Month  data  -  The  calculation  of  the  Root 
Mean  Square  error  will  be  discussed  in  III.  E.  (Monterey  Station.) 

•  Forecast  the  tide  for  the  next  24h  using  “t_predic”  and  the  structure  of 
constituents  that  resulted  from  the  detiding; 

•  Low  pass  fdter  the  residual  current  to  1 8h; 

•  Average  the  residual  current; 

•  Take  the  last  24h  wind  readings; 

•  For  each  hour  in  the  future  h,  get  the  forecast  wind  from  that  hour  until 
t=0  and  the  wind  readings  before  that  until  h-24; 

•  Low  pass  fdter  both  to  1 8h; 

•  Do  the  average  for  each; 

•  Take  the  difference  of  the  winds; 

•  Use  the  transfer  function  to  get  the  current  associated  with  the  wind 
difference; 

•  Sum  the  residual  current,  the  tide  and  the  wind  current. 

Figure  28  is  an  example  of  the  obtained  result,  in  a  day  where  a  major  wind 
change  took  place  and  the  wind  influence  is  noted.  The  resolution  is  hourly. 
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(m/s)  (m/s) 


01  /Jan  u  rms=0. 060447 


01 /Jan  v  rms=0.047992 


Figure  28.  One  day  forecast  with  tide  treatment  -  The  wind  changing  its  influence. 
(Monterey  Station.) 
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2. 


With  Persistence 


The  procedure  sequence  is  the  same  as  the  above  but  with  the  major  difference  of 
not  removing  the  tidal  influence.  This  method  uses  the  last  24h  current  readings  as  the 
current  for  the  next  24,  summed  with  the  wind  difference  driven  current.  Example  results 
using  the  persistence  model  are  given  in  Figure  29. 


01  /Jan  u  rms=0. 057852 


01 /Jan  v  rms=0.06595 


Figure  29.  One  day  forecast  with  persistence  -  The  wind  changing  its  influence. 
(Monterey  station.) 
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D.  TUNING 


The  wind  interval  to  account  for  the  difference  between  forecast  and  real  is 
chosen  after  several  different  runs  bases  on  shifting  the  future  wind  interval  for  several 
periods  forward.  The  one  with  least  error  is  the  one  designated,  where  the  wind  interval 
ends  at  tp  (Figure  30).  Several  values  around  the  transfer  function  were  also  attempted. 
The  best  results  are  obtained  with  a  50°  clockwise  rotation,  and  with  a  magnitude 

Cwiltd  =0.01  SU5  +0.018.  (12) 


_ 

wind  to  take  the  difference 

last  24h  wind 

l 

t 

t 

1 

to 


Figure  30.  The  difference  between  winds  -  Several  positions  of  the  future  wind  interval 
were  attempted  for  a  given  time  of  prediction  (tp). 

E.  COMPARISON 

RMS  is  calculated  on  a  daily  basis  with  and  without  introducing  the  wind  .  RMS 
is  based  on  the  difference  between  our  calculated  surface  current  for  each  hour  and  the 
actual  CODAR  reading  at  that  time.  This  means  that  we  are  considering  CODAR  as 
“true.”  We  must  take  into  account  the  CODAR  intrinsic  variability  when  analyzing 
results.  The  end  RMS  result  for  each  component  u  and  v  is: 


with  the  total 


X  ( reading  -  forecast )2 


RMS 


component 


N- 1 


(13) 


RMSlolal  =  pMSf  +  RMSf  .  (14) 

The  results  are  given  in  Table  11,  considering  a  normal  distribution. 
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(m/s) 

average 

Variance 

Standard 

deviation 

95%  conf 

(+/-) 

99%  conf 

(+/-) 

Id 

> 

RMS 

0.1593 

0.0022 

0.0471 

0.0923 

0.1213 

o 

£ 

<L> 

J-H 

0) 

With  wind 

RMS 

0.1667 

0.0028 

0.0529 

0.1036 

0.1363 

H 

No  wind 

0) 

o 

RMS 

0.1427 

0.0023 

0.0482 

0.0945 

0.1242 

£ 

0) 

With  wind 

c/3 

u 

0) 

RMS 

0.1495 

0.0028 

0.0528 

0.1035 

0.1360 

Oh 

No  wind 

Table  11.  Root  Mean  Square  comparison  -  Average  daily  values.  (Monterey  station, 
Jul  06  to  Dec  07.) 


The  statistical  results  for  the  18-month  study  period  show  that  the  wind  influence 
is  one  order  of  magnitude  inferior  to  the  total  error.  Use  of  persistence  or  tide  treatment 
yields  practically  similar  results,  with  simple  persistence  of  today’s  observed  low  pass- 
filtered  currents  giving  slightly  better  results. 

Looking  at  other  station  results,  station  MO  gets  the  least  error  results,  followed 
by  Ml,  probably  due  to  the  repetitive  nature  of  the  flow  in  that  position,  or  perhaps 
because  it  is  more  protected  from  unaccounted  variations.  Results  distributed  by  station 
are  shown  in  Figure  31.  There  is  a  visible  relation  between  error  and  stations,  showing 
that  the  error  causes  are  probably  associated  with  the  position  and  its  typical  patterns  of 
wind,  currents,  the  radar  station  itself,  and  other  features.  This  will  be  studied  more 
closely  in  Section  III.F.  Interestingly,  the  distributed  results  show  clear,  positive 
correlation  between  forecast  error  and  the  mean  station  current  speed  or  current 
variability,  which  is  consistent  with  the  previous  findings  of  Garfield  et  ah,  2007. 
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Figure  3 1 .  Relation  of  RMS  error  and  the  stations  -  Top  is  magnitude  of  the  current  vs. 

error.  Bottom  is  variability  vs.  error.  Based  on  the  tide  plus  mean  with 
transfer  function  model  forecasts. 
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Seasonality  does  not  seem  to  play  a  major  role,  but  the  small  quantity  of  data  does 
not  allow  a  definitive  conclusion.  The  four-season  results  for  the  Monterey  station  are 
shown  in  Figure  32. 
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Figure  32.  Relation  of  RMS  error  and  seasons  -  Top  is  magnitude  of  the  current  vs. 

error.  Bottom  is  variability  vs.  error.  (For  Monterey  station,  based  on  the 
tide  plus  mean  wind  transfer  function  forecasts.) 

The  long  records  from  the  Monterey  station  can  be  used  to  illustrate  the  temporal 
behavior  of  the  forecast  error  and  the  influence  of  the  wind-based  corrections.  The  daily 
RMS  current  speed  errors  are  plotted  in  Figure  33  with  and  without  the  wind-based 
corrections.  Almost  everywhere,  the  wind- treated  forecasts,  blue  and  black,  are  lower. 
Red  and  green  tips  are  easily  seen  constantly  at  the  top,  despite  the  small  difference. 
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RMS  compare  -  Wind  influence 
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Figure  33.  Zoom  of  a  plot  with  the  different  calculated  RMS  -  Wind  influence,  “p”  is 
persistence,  “nw”  is  no  wind.  (Monterey  station.) 

Zooming  in  on  portions  of  the  record  (Figure  34)  allows  us  to  see  the  nature  of 
some  types  of  errors.  As  one  example,  sometimes,  not  often,  the  RMS  associated  with  the 
tides  treatment  forecasts  goes  up  with  no  apparent  cause,  probably  the  result  of  a  bad 
frequency  fitting. 
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Figure  34.  Zoom  of  a  plot  with  the  different  calculated  RMS  -  Possible  tide  problem, 
“p”  is  persistence,  “nw”  is  no  wind.  (Monterey  station.) 

F.  ERROR  CAUSES 

If  we  take  into  account  the  average  current  velocities  in  this  area,  0.17m/s,  errors 
of  0.15m/s  represent  almost  the  same  magnitude.  Nevertheless,  considering  the  main 
objectives  of  the  application,  having  a  measure  of  the  uncertainty  is  a  good  result. 

If  we  have  a  distress  in  a  position  where  the  rescue  mean  takes  T  time  to  fulfill  the 
transit,  with  a  C  uncertainty  current,  after  the  transit  the  search  area  will  be  at  least  a 
circle  of  radius  TxC  =  M  meters.  Reducing  the  C  will  proportionally  reduce  M,  which 
means  that  it  will  take  less  time  to  find  the  target. 

The  smaller  the  error,  the  bigger  the  odds  of  saving  a  life  or  reducing  an  oil  spill’s 
impact.  With  this  in  mind,  continuation  of  this  research  should  be  a  priority  of  the  naval 
sciences. 
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1.  Accuracy  Expectancy 

In  order  for  us  to  have  an  idea  of  what  to  expect  in  terms  of  accuracy,  we  need  to 
layer  known  uncertainties  related  to  the  surface  currents. 

a.  CODAR  (From  Wright,  2008) 

In  a  good  scenario,  we  will  have  approximately  0.08m/s  accuracy. 

b.  Winds 

The  influence  is  unknown,  but  considered  to  be  small,  according  to  our 
transfer  function  result.  The  reading  error  multiplied  by  one  order  of  magnitude  below  the 
induced  current  error. 


c.  Tides  (Pawlowicz  et  ah,  2002) 

Using  the  error  output  of  t_tide  and  averaging  the  errors  of  the 
constituents,  for  one  month,  we  have  an  uncertainty  of  around  0.020m/s.  If  we  use  3  days 
the  uncertainty  is  almost  the  same,  0.022m/s. 

d.  Growing  Surface  Waves  (Ardhuin  et  ah,  2003) 

This  phenomenon  may  “induce  stresses  opposed  to  the  wind  stress  for 
wave  growth  stages  that  may  represent  up  to  10%  of  the  wind  stress  for  short  fetches.” 
For  a  typical  current,  and  according  our  transfer  function,  this  represents  an  uncertainty 
of  about  0.02m/s. 


e.  Eddy  Surface  Velocity 

In  this  study,  a  typical  eddy  speed  variation  over  24h  is  considered  to  be  at 
least  0.02m/s.  Eddies  are  highly  variable  but  usually  typical  of  certain  locations.  This 
uncertainty  depends  on  the  physical  characteristics  and  if  the  eddy  is  influencing 
movement  near  the  core  or  the  out  skirts.  Observing  the  relative  path  is  also  important. 

f  Others 

Other  kinds  of  possible  interferences  were  considered,  but  not  included: 
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(1)  Stokes  Drift.  For  gravity  waves: 

Us  =  a2cok  (15) 

where  a  is  the  amplitude,  to  the  frequency  and  k  the  wave  number  (Radko,  2007,  Ocean 
Dynamics  class  notes  NPS). 

This  yields  for  a  lm  amplitude,  10s  period,  10m  length  wave  an 
induced  velocity  of  0.39  m/s,  which  means  an  order  of  magnitude  similar  to  the  total 
RMS  error.  Studies  show  that  it  may  be  opposed  and  cancelled  by  a  Eulerian  flow  (Xu 
and  Bowen,  1993;  Monismith  et  ah,  2006).  Works  are  still  in  progress  to  evaluate  the 
total  waves  i-duced  flow. 

(2)  Inertial  Movements.  These  play  a  role  but  depend  on  the 
local  flow  characteristics.  They  do  not  seem  to  play  a  major  influence  in  the  Monterey 
buoy  area  circulation.  Its  surface  currents  may  be  over  0.1  m/s,  again  the  order  of  our 
RMS  error.  We  have  to  take  into  account  that  if  using  yesterday’s  current  for  tomorrow, 
this  error  will  be  minimized.  If  we  use  tide  treatment  in  a  zone  of  important  inertial 
motions,  care  should  be  taken  before  averaging  for  the  mean  flow. 

(3)  Other  Non  Linearities.  (Eddies,  larger  periodicity  waves, 
etc.)  All  of  these  influences  exist  and  are  a  part  of  the  flow.  The  use  of  persistency 
already  takes  into  account  enough  of  its  magnitude. 

The  total  sum  of  the  uncertainties  is  0.14m/s.  Once  these  values 
were  quantified  in  a  conservative  way,  this  means  that  the  predictable  accuracy  is  not 
expected  to  be  better  than  the  order  of  our  actual  results. 

2.  Abrupt  Changes 

The  plotting  of  RMS  error  against  significant  wave  height  (SWH)  and  current 
modulus  (following  Garfield  et  ah,  2007)  is  almost  intuitive  and,  by  visual  inspection, 
points  out  possible  error  causes.  Wave,  wind  and  current  variations  are  shown  as  time 
series  in  Figure  35,  together  with  the  RMS  forecast  current  error  time  series.  Whenever 
an  abrupt  change  exists  in  the  current  speed,  RMS  seems  to  go  up.  The  wind  acts  against 
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this  rise,  decreasing  the  RMS  each  time  it  couples  the  change  in  the  current.  Wave  height 
seems  not  to  have  a  significant  consistent  impact.  This  impression  is  clearly  supported  in 
the  zoomed  time  series  presented  in  Figure  36. 


RMS  forecast  with  persistence  and  wind 


Figure  35.  Full  RMS  plot  versus  SWH,  current  and  wind.  (Monterey  station,  Jul  06  to 
Dec  07.) 


60 


RMS  forecast  with  tides  treatment  and  wind 


RMS  forecast  with  persistence  and  wind 


Figure  36.  Zoom  of  RMS  versus  SWH,  current  and  wind  -  High  error  seems  associated 
with  large  current  variation.  Top  with  tides  treatment.  Bottom  with 
persistence.  (Monterey  station.) 
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IV.  CONCLUSIONS 


This  study  analyzed  long  records  of  HF  radar-derived  surface  currents  and 
mooring-based  wind  observations  for  the  region  offshore  Monterey  Bay,  California.  The 
growth  of  HF  radar  instrumentation  and  the  resulting  surface  current  mapping  data  has 
enabled  a  number  of  scientific  and  public  safety  capabilities.  In  many  instances, 
particularly  with  respect  to  search  and  rescue  or  spill  mitigation  applications,  surface 
current  forecasts  are  critical  in  addition  to  real  time  observations.  To  make  accurate 
forecasts  requires  knowledge  of  background  flow  patterns  plus  variable  flow  patterns 
related  to  tides  and  wind-driven  currents.  The  focus  here  has  been  on  created  and  testing 
models  for  wind-driven  surface  currents  that  are  specific  to  the  near-surface  (~0.5m) 
depths  of  the  HF  radar  observations. 

The  Monterey  Bay  region  contains  a  wealth  of  surface  observations  from  both  HF 
radar  and  meteorological  buoys.  One  site,  in  particular,  surrounding  the  NOAA/NDBC 
buoy  46042  location  has  provided  an  eighteen-month,  continuous  set  of  currents  and 
winds.  Those  data  were  used  to  develop  an  empirical  transfer  function  between  wind  and 
HF  radar-derived  current.  Shorter  data  sets  from  other  buoy  sites  both  inside  and  outside 
Monterey  Bay  were  used  to  confirm  results  from  the  buoy  46042  location.  In  addition  to 
the  forward  problem  of  creating  a  best-fit  transfer  function  from  wind  to  current,  the 
recommended  parameters  were  adjusted  slightly  by  determining  which  values  provided 
the  most  accurate  24-hour  forecast  of  HF  radar-derived  surface  current  at  the  buoy  46042 
location.  The  results  of  these  wind  versus  surface  current  analyses  from  the  1.5  year 
dataset  from  five  different  stations  lead  to  the  following  conclusions  and  requirements. 

A.  THE  WIND  INFLUENCE 

The  wind  induced  current  is  on  the  order  of  2%  of  the  total  wind  speed,  which  is 
consistent  with  historical  values.  Despite  this  magnitude,  the  influence  of  introducing  an 
explicit  wind-driven  current  correction  to  the  existing  mean  plus  tidal  current  forecast 
models  leads  to  around  10%  of  increased  accuracy  in  the  forecasts.  The  best  results  are 
obtained  with  a  50°  clockwise  rotation,  and  with  a  magnitude  Cwind  =  0.0 18U5 +0.018 . 
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Although  this  study  has  shown  a  clear  benefit  to  including  a  wind-driven  correction,  this 
correction  is  still  below  the  level  of  the  forecast  uncertainties  that  derives  from  other 
uncorrected  factors. 

B.  RESULT  CAPABILITIES 

The  overall  result,  including  the  wind-driven  current  correction,  shows  that  HF 
radar-derived  surface  current  mapping  data  are  able  to  forecast  surface  currents  in  the 
next  24h  with  errors  of  about  0.15  +/-  0.1  m/s  with  95%  confidence.  When  significant 
changes  in  current  happen  without  the  wind  influence,  approximately  one  day  of  high 
RMS  forecast  error  is  expected.  The  largest  and  most  persistent  error  events  for  the  test 
period  from  July  06  to  December  07  led  to  a  four-day  RMS  error  over  0.2  m/s. 

C.  NEXT  STEPS 

The  goal  here  has  been  to  develop  improvements  to  the  short-term  surface  current 
forecast  model,  with  primary  focus  on  included  the  effects  of  predicted  wind  changes. 
Several  other  factors  that  are  expected  to  influence  the  accuracy  of  24-hour  surface 
current  forecasts  were  pointed  out  in  Chapter  III.  The  results  here  show  that  it  will  be 
necessary  to  understand  and  incorporate  more  of  those  influences  if  improvements 
beyond  the  10%  level  achieved  here  are  desired.  The  qualitative  assessment  of  error 
sources  presented  in  the  last  chapter,  combined  with  previous  work  on  HF  radar  surface 
current  validation,  suggests  that  errors  in  the  measurements  themselves  may  still  be  the 
biggest  limitation  to  improved  forecasts. 

In  addition  to  work  that  is  ongoing  to  improve  HF  radar  surface  currents  and  error 

tracking,  this  study  points  to  the  need  to  begin  to  incorporate  the  most  important  aspect  of 

HF  radar  observations  into  the  forecast  process.  That  is,  spatial  infonnation  from  the 

surface  current  mapping  data  should  be  incorporated.  At  this  time,  all  information  about 

the  sub  tidal,  tidal,  and  wind-driven  surface  currents  that  goes  into  the  forecast  models 

derives  from  time  series  information  from  a  single  location.  HF  radar,  on  the  other  hand, 

provides  a  map  of  surface  currents  every  hour.  Observations  show  those  maps  contain  a 

good  deal  of  horizontal  variation  in  the  surface  currents,  which  is  the  fundamental  reason 

that  a  network  of  HF  radar  systems  is  necessary  for  the  transport  applications  supported 
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here  compared  with,  for  example,  a  single,  moored  current  meter.  The  next  generation 
forecast  models  should  account  for  advective  effects.  The  use  of  spatial  integration 
between  one  calculated  point  and  its  immediate  neighbors  may  allow  gross  errors  to 
emerge  more  clearly. 
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